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THE METEOROID EWIROEMENT AND STRUCTURAL RELIABILITY 

1,0 jJTRODUCT-lOH 


Primarily as a result of the weight limitations on Apollo, a com- 
prehensive survey of the majority of all. the- articles ever- published, on 
meteors, "by competent authorities in the field, was undertaken in an 
attempt to establish the validity -of the 1957 Whipple environmental - 
model. This period of review culminated in a visit to the Harvard College 
Observatory and the Smithsonian Institution Astrophys ical Observatory for 
consultations with Drs. Whipple,. McCrosky, Cook, Jacchia, Hawkins, and 
Southworth, 

As the result of all this activity, it was possible to redefine- the 
meteoroid environment based on’ fact tempered with the respected opinions 
of the above named experts. The .new. environment, defined in this paper, 
was compiled by the Apollo management and 'support groups 'as an interim 
revision. Further revision is inevitable as the data from current and 
proposed meteoroid sampling programs come to hand. 

The position regarding penetration theories has hot altered much 
in the interim between the Apollo study phases and the .present time. 

There is still no concrete proof that Bjork’s theoretically deduced • 
equation is any more valid than the empirical Summers -Charters equation 
at velocities within the meteoric range. The promise of really high- . 
speed impact devices, with controllable masses and velocities, has not 
materialized to date;, hence the solution of this .vexing problem has been 
delayed. ■ 

This paper attempts to lay the foundations of a clearer understanding 
of the problem facing those who have to design for a seemingly nebulous 
environment. Only those factors truly important in this regard have been 
included in the appropriate sections dealing with the whole story. 

The appendix contains worked examples on the evaluation and design 
of structure for both the no penetration and the controlled penetration 
cases. Material selection is dealt with purely from the standpoint of 
the two penetration theories considered in this paper. 

It has not been possible to include any data specifically pertinent 
to nonmetallic materials under’ ’impact. Current practice Has been to 
assume that the results obtained for metallic targets are applicable. 
Another omission is the subject of micro -meteoroid erosion and it is 
hoped that a future paper will rectify these shortcomings. 
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SYMBOLS 


Area, sq. meter or sq. ft. 

A constant 
Confidence level 

Sonic velocity in target material, km per sec 

.Diameter 'of projectile, cm 
QJrue failure rate 
Most probable failure rate 
Luminous flux, vis vial ergs per sec 

Penetration resistance efficiency factor; a constant 
F 

Ratio — ~ 

F 

Luminosity coefficient,, sec per km 

Luminous efficiency 

Luminosity, foot candle 
Meteor magnitude 
Meteor visual magnitude 

Meteor mass, grams 

Meteor mass at zero visual magnitude, grams 

Mass of projectile, grams 

Meteor flux, number per sq. meter day 

Meteor flux at' a given magnitude, number per sq. meter day 
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Meteor flux at zero visual magnitude, number per sq. meter day 

Np Penetrating meteoroid flux, number per sq. meter day 

W Total meteoroid flux, number per sq. meter day 

n Velocity index; number of trials 

"til. 

P^ Probability of r event 

p Probability of single event in a single trial; penetration, 

depth, cm 

1 (1 - p) 

R True reliability 

R Most probable reliability 

t Duration of meteor trail, secs 

T Equivalent single sheet thickness, cm 

t „ Thickness of individual sheet, cm 

f 

V Meteor velocity, km per sec 

V Meteor velocity at zero visual magnitude, km per sec 

o 

V Relative impact velocity of projectile, km per sec 
P 

V Relative impact velocity of meteoroid, km per sec 

a Meteoroid flux-mass constant, gm impacts per sq. meter day 

0 Meteoroid flux -mass index 

T| Earth shielding factor 

3 

Meteoroid density, gm per cm 
p ' Projectile density, gm per cm^ 



Target material density, gm. per cm" 


T Mission time in environment, days 


2.0 METEOR THEORY 


A precise definition of all the factors involved in the formula- ' 
tion of a design meteoroid environment is not currently available for 
the entire particle spectrum. While it is true that steps are being 
taken to obtain the data that has been lacking, the present generation 
of space-vehicles will he designed long before the results are available. 
Fortunately, it is possible to synthesize a model environment based on 
the known facts, using reasonable assumptions whenever necessary. An 
outline of the fundamentals involved in this task is presented in the 
ensuing text, without recourse to the history of meteor astronomy which 
may be obtained from references 1 and 2. 


2.1 METEOR MAGNITUDE ' 


The apparent magnitude of a star or a meteor is a measure of its 
luminosity as observed from the earth, relative to a reference stellar 
zero magnitude. A logarithmic comparative scale was adopted in 185 O 
based on the equivalence of a magnitude interval of unity to a luminosity 
ratio of 2. 512. That is 



2.512 



( 2 . 1 ) 


where and are the magnitudes and and the correspond- 

ing luminosities. 


In logarithmic form, equation ( 2 . l) is 


log 


10 L„ 


°- 4 ( M a -*0 


( 2 . 2 ) 


The scale is inverted in that increasing positive magnitudes are 
indicative of fainter objects. 
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A "visual" magnitude scale is most commonly used in metepr studies 
in spite of the fact that photography and radio astronomy have long 
since replaced the visual and optical observational methods. The more 
recent determinations of meteor flux were obtained in terms of photo- - 
graphic and radar magnitudes which were then converted into visual 
magnitudes by means of predetermined correlation constants.,, " 

Whichever scale is used, it functions as a means of classification 
of the meteoroid population in terms of flux, velocity, and, more 
recently, density. 

A relationship between the absolute visual magnitude of a meteor, 
M v , which is defined as . the magnitude at a distance of 100 km in the 
direction of the zenith and the luminous flux, X, in visual ergs per 
second, has been derived by Opik, reference 4 , 

This relationship is, 


log l0 I « 9-72 - 0 . 4 M v (2:5’) 

based on a stellar magnitude of -26.72 for the sun and the energy 
distribution of solar radiation. 


2.2 THE RELATIONSHIP BETWEEN METEOR FLUX AND VISUAL MAGNITUDE 

The number of meteors dN^ appearing in the atmosphere between 
visual magnitudes and is given by the incremental law. 


dN _ = a. r Mv . dM 
M v 


( 2 . 4 ) 


where a is a constant and r is the ratio of the increase in the 
number of meteors between magnitudes M and + l. 

Integrating equation ( 2 . 4 ) between the limits «= -« and M 
gives the cumulative law.. In logarithmic form, this is, 

. lo ^10 h m = M V 1Og 10 r + f og 10 a - l0g l© ( ln r >) 


- M v • loE 10 r + log 10 H o 


( 2 - 5 ) 
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where BT is the number of meteors with M = 0 and "brighter, and -r 
o v 

is the same as in the incremental law. 

The cumulative law is the "basic concept in "the. definition of the 
meteoroid environment and there is no disagreement amongst astronomers 
as to its validity. 

Values of r vary "between 2. 0 and k. 0, and of II , "between 

6 *5 ^ 

1.0 X 10 and 1.9 X 10 , as a result of the different observational 

techniques used, the uncertainties in the observations, and in the 

correction factors involved. 


McKinley, in reference 2, gives the flux -magnitude relationships 
favored "by Watson, Millman, and Hawkins and Upton, as: 


a. 

Watson (1958). - 

. lo S 10 

= 5. 88 + 0. 40 M 
M V 

for 

-3 < M <10 
■— v — 

b. 

Millman. - 

log 10 

= 6/0 + 0.57 M 
. M v 

for 

-10 < M < 0 

_ y — 



log io 

JSL. = 6. 0 + 0. 50 M 
M V 

for 

tr\ 

Vi 

VI 

0 


- 

^■ Og io 

U « 6. 5 + 0. 40 M 
M ■ V 

for 

IA 

<1* 

1 A 

H 

O 

c. 

Hawkins and Upton. - 

1O ®10 

\ = 5.27 + 0.537 \ 

for 

-! < < 10 


The limiting magnitude of +10 in all of these expressions is also 
the limit of the reliable radar, observations to date and until systematic 
observations can be extended beyond this limit, it is necessary to assume 
that they are applicable up to M = + 15 . 

One is forced to agree with "Hawkins and Upton's contention that 
as -> a. , r must approach a value between 1 and 0 , to avoid the 

inference that there is an infinite number of meteoroids in the universe. 

A plot of as a function of M , reproduced from reference 2, 

is presented in figure 1 to show the relative significance of the three 
relationships listed above on the Apollo spacecraft. The number of 
meteors of zero magnitude and brighter is based on the mean hourly flux 
rate of the sporadic meteor activity alone. This mean hourly rate is 
obtained from the results of several years of photographic observations, 
-and for the Super-Schmidt camera collection area of 5 .>980 km 2 at an 
altitude of 90 km. 

* K 
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Meteor flux may be quoted as the number per square meter per hour 
or per day, or as the number over the entire .surface area of the- earth 
per day. In the latter case, the surface area of the earth is calcu- 
lated at an altitude of 100 km. 

2, 3 THE CCffiRELAT-lOE OF METEOR MAES MB FLUX 

Due to the evidence that most meteors break up upon entry into the 
earth’s atmosphere, the most reliable method of calculating the mass of 
a meteor is by the relationship between the luminosity and the kinetic 
energy of a meteor, 

1. 1 i k^.mV 2 (2.6) 


In this expression, I is the luminous flux which is related to 
the meteor magnitude M by equation (2..3)> t - is the time' during whiel: 

the meteor is visible, k is the luminous efficiency factor, and — mV^ 

is the kinetic energy of 'the-ipeteor. 

The value of quoted is generally based on the photographic- 

magnitude corresponding to the point of maximum luminous flux of the 
meteor trail. As a result, equation (2.6) may only be used for meteors 
which exhibit very little change in magnitude over the length of the 
trail. 

Most meteors exhibit a variation in luminous flux along the length 
of the trail, and the time -histories of three types are shown in figure 

In the general case, equation (2.6) must be replaced by, 

00 

I «dt = 
t 



m 


i V 2 


(2.7) 


where the integration limits refer to the luminosity between the time 
when the mass is m and + «, when the meteoroid enters the atmosphere. 


The terminal mass is generally zero for most photographic and • 
radio meteors, hence it is possible to calculate the original mass of 
the meteoroid, m , by integrating I- dt between the limits + and 
that is , 



I - dt 



m 

00 


V 2 


( 2 . 8 ) 
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' Until "very recently, it has not been possible to determine the 
luminous efficiency directly for various reasons. Most of the determi- 
nations of meteor mass have relied upon theoretical considerations of 
the lavs of atomic physics, and estimates of have ranged between 

2 k 

2 X 10 and 2 X 10 , Foremost in this field has been the work of Opik, 

k l 

and in reference 4 , he has published a table of values of k ® — for 

o 

various meteor velocities ranging from, 5' to 84 km/sec. 


This- quantity, k Q , is sometimes referred to as the luminosity 


coeffic ient. 


The dependence of k on meteor velocity to some power n is 

recognized amongst astronomers and it is thought that for the brighter 
photographic meteors, say up to magnitude * +4. 0, n- is-., unity; A ■ 
relationship for the fainter photographic and radio meteors remains to 
be established. 

In the calculation of meteor mass, the astronomers actually use 

k rather than k, so that equation (2.8) is rewritten as, 
o 1 

p-bo k , 

/ I * dt = “ • (2., 9) 

M-oo 


A zero magnitude mass is established using the appropriate velocity, 
and then the mass value for increasing or decreasing positive magnitudes 
is obtained from the assumption of a mass ratio of 2.512 per unit 
magnitude step. It is also assumed that the mass of a meteor is directly 
proportional to its luminous flux so that increasing positive magnitudes 
correspond to decreasing masses. 


The expression linking mass and .magnitude is. 


from which 


m = 


10 


K 
o* 4 m 


v 




m 2 


= 10 


- m \ (m , - m \ 

v =2 = 2 .- 512^ V=1 T=2 / 


( 2 . 10 ) 


A summary of the zero magnitude mass, velocity assumed, luminosity 
coefficient used when known and meteoroid density is as follows: 
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Watson (194-1) — 


ni Q = 0.25 gm. 


V * = 55 km/sec 


• k o 3 ' 


P m = 3 . k 


m o = 1. 25' gm 


■V =4-0 km/sec 


k = 8. 5 x IQ”- 5 see/km p =4 gm/cc 
o m 


= 25« 0 gm 


V a 28 km/sec 
o ' 


k * 1.24- x 10"^ sec/km p =0.05 gm/cc 


Whipple & Hawkins — m = 0. 15 gm 
(1958) 


= 50 km/sec 


k = 3. 5^- X 10 sec/km p = ? 
o ' m 


McCrosky &. Posen — m = 0. 17 gm 

(1961V 


V q - 28,9 km/sec 


k = 8.5 X 10 -5 sec/km P = ? 

r\ * m * 


m o = 1.29 gm 


V = 4-2 km/sec 
o ' 


k = 1. 06 X 10“^ sec /km p = < 0. 1 gm/ee 


m = 0. 29 gm 
o 


V q = 42 km/sec 


k a 4. 98 X 10 sec/km- P =3-4 gm/ce 


0.1 Whipple (1961) - 


m Q = 2. 5 gm 


V =28 km/sec 
o ' 


k = 1,24 x 10” sec/km P =0.5 gm/cc 
o m Q ' 
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McCrosky & Cook — m = 2.0 gm V = 28 km/sec 

(19^7 

k = 2.29 X 10 sec /km = 0. 6 gm/ec 

The immensity of the problem of assigning the zero magnitude mass 
is exposed by this historical resume of the varying estimates’ of the 
astronomers. The acquisition of data from the Trailblazer projects of 
1962 ^ when artificial meteors of known mass were fired into the atmos- 
phere and photographically observed from ground stations, has resulted 

in a directly calculated value of k Q - 2.29 X 10~^ sec/km as determined 

by Cook, Jacchia, and McCrosky. 

This value results in the zero magnitude mass calculated by McCrosky 
and Cook, ref. 6, shown in the last line of the above summary. It is \ ' 

interesting to note that the value of k Q obtained by Hawkins and Upton, 

using an ionization efficiency, the counterpart of luminous efficiency 
in radio astronomy, is no more than one order of magnitude less tharwfehe 
Cook, Jacchia, and McCrosky value mentioned previously. 

The relationship between meteor mass and the particle flux is obtained 
from equations (2.5) and (2.10) and is of the form, 

F = a-m* 3 (2.11) 

where a and p are constants, N is the impact rate and m is the 
meteoroid mass. 

When N is measured in impact rate per square meter per second and 
m is in grains, the values of a, and P are as follows: ' 



Source 

a 


P 

Mass Range (gm) 

a. 

Meteors 

■ 5- 715 X 

10- 1 ?- 

-1.55 

10" 5 to 10° 

b. 

1957 Whipple 

1. 518 x 

10 " 12 

-1.00 

10' " 9 to io° 

c. 

1961 Whipple 

5.012 X 

10~ 15 

-1. 186 

10" 9 to 10° 

d. 

1956 Watson 

5- 129 x 

10” 15 

-1. 00 

10 “ 9 to 10° 

e. 

0. 1 Whipple 1957 

1.518 X 

10 “ 15 

-1. 00 

10" 9 to 10° 

f,. 

0.1 Whipple 1961 

5-012 X 

1Q- Ik 

-1. 186 

I0~ 9 to 10° 


The values for (a), (b), and (c) were obtained from reference J , (d) is 
based on Watson's value of m Q = .25 gm, the flux -magnitude relationship 
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of section, .2. 2, .and 9° X' 10 ' meteors' per 'day at = +5. The zero 

magnitude mass of (t>) and (c) is 25 gins, whereas (e) and '(f), acting on 

• 1 ■ . ■' 

Whipple's- suggestion in- reference 7> employ -Jq of this' value. • 

Log -log plots of some of these flux -mass relationships are given 
in figure 3- ’ 

Figure 3 also indicates the. change in slope which results if the 
satellite data .of reference 8 is representative of the environment in 
the entire earth-moon system. . . , 

Whipple's contention of a concentration- of dust around the earth 
which follows an approximately, inverse 1. 4 power law in distance from 

the earth's, surface to about lCr km, reference 7, is shown, in figures 4 
and 5* 

2.4 .METEOR VELOCTFT 


One of the facets of the meteoroid environment that has received a 
great deal 'of attention from the astronomers is the problem of ‘ the 
limiting maximum velocity. Since the beginnings of the scientific obser 
vat ions the origin of meteors has been the subject of controversy. 
Basically, if meteors are members of the solar system, they must have 
closed orbits and the limiting heliocentric velocity at the earth's 
distance from the sun for this case is 42 km/sec. 

As the earth's heliocentric velocity is approximately 30 km/sec, 
in the extreme "case of a retrograde meteor in the plane of the ecliptic, 
the observed velocity relative to the ■'earth cannot exceed 72' km/sec. 
Whipple in reference- 6 gives a value of 73 km/sec" including a slight 
additional effect due to the earth's gravitational potential. 

Primarily as the -result of radio velocity- observations, -astronomers 
are now- in agreement that practically all the meteors which, have been 
observed to date are members -of the solar system. McKinley, In ref erencs 
2, states that of 11,000' radio meteors, only '32 had a velocity in the 
range of 75 to 79 km/sec and there were none exceeding 80 km/see. He 
concludes that some of these' could he explained by the 5 percent error 
in the reduction process. 

Both Lovell reference 1, and Jacchia and Whipple, reference 5? in 
their radio arid r photographic velocity determinations of several hundred 
meteor sightings could find no more than 1 percent that had velocities; 
exceeding the critical value. 
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Tims it seems conclusive that the upper limit of the geocentric 
velocity is approximately 72 km/sec. 

In the case of a direct meteor in the ecliptic, the observed velocity 
relative to the earth must, - hy the same token, be approximately 12 km/sec. 

Within this range of 12 to 72 km/sec for the photographic and radio 
meteors, the average velocity has been calculated by various investigators. 
Thus Whipple, .in reference 1 j , gives an average of 28 km/sec for photo- 
graphic meteors. In reference 2, McKinley reports that the average ob- 
served velocity obtained for the 11,000 radio meteors is 44,7 km/sec, while 
the radio results for two sets of nearly 2,000 radio observations made 
at Jodrell Bank in England were 44.7 and 59-2 km/see. For the photo- 
graphic meteor observations, McKinley reports that the average entry 
velocity of 621 meteors reduced by Jacchia, Hawkins, and Southworth is 
55-6 km/sec, and 3^-3 km/sec for the 2,433 velocities determined by . 
McCrosky. 

The characteristic velocity distribution plots for radio and 
photographic meteors are shown in figures 6 and 7- 

A program to determine the average entry velocity using McCrosky; 's 
data for the 2,433 photographic meteors, as well as the average per 
magnitude interval, has recently been concluded within the Spacecraft 
Technology Division, reference 10. 

The preliminary results are shown in figure 8, the average in this 
case being 32.7 km/sec. 

Velocity variation in terms of visual magnitude has received scant 
attention, and Whipple's attempt to establish a scale, reference 3, 
remains the only one published by the experts in this field. 

The entry velocity at the edge of the earth's atmosphere cannot he 
less than 11 km/sec for a particle starting from rest within the solar 
system. Whipple, in reference 3> says that higher velocities can occur 
for nearly circular orbits by inclination of the orbital plane to that 
of the earth. For this reason, he .arbitrarily chose a mean value of 
15 km/sec for the smallest particles and then .assumed a linear varia- 
tion between this value and the average velocity for the photographic 
meteors. 

Figure 9 is a plot of geocentric velocity as a function of visual 
magnitude for Whipple's estimate, together with the author's estimates 
based on the average velocities cited in this section. 
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2. 5 THE density of meteoroids 


The examination of meteorites found on the earth, which are meteor- 
oids that have survived the passage through -the earth's atmosphere, has 
established the fact .that! 

a. Their chemical composition is unlike any rocks found on the 
•earth. 

b. About 90 percent are stony in character, with an average 
density of 3*5 gm/cc, 

c. About 10 percent are iron-nickel with an average density of 
7, 8 gm/ec. 

The 'identification of meteor orbits which do not exhibit any great degree 
of eccentricity and, in addition, have aphelion points - inside the- orbit 
of Jupiter, has led to the belief that meteorites originate . in the . 
asteroidal belt between Mars and Jupiter. : 

Meteorites are considered as part of the total meteor flux ahd are 
commonly called asteroidal meteors. The percentage of the total- flux 
that is asteroidal varies between 10 percent, as determined by Jacchia 
and Whipple in reference based on the orbital characteristics of 
meteors, and no more than 2 percent, based on the photographic evidence 
of fragmentation. 

This latter figure is the opinion of both McCroSky and Cook of the 
Harvard Meteor Department, and it is considered applicable within the 
range of visual magnitudes 0 to +5. 

Radio meteor studies to magnitude +10 currently indicate that 
fragmentation is again predominant so that by the. same reasoning it 
seems likely that the asteroidal content is Very small. 

Considering the evidence of meteorites, and fireballs, it is a 
fair assumption that the asteroidal content of the total flux increases 
for magnitudes- brighter than Zero. • 

1 

The density of the asteroidal component of the total flux, because 
of the preponderance of stony meteorites, is generally given as 3*5 gm/cc 

Association of shower meteor orbits with comet orbits was first 
noticed in 1862, and since that time 12 meteor streams have been, found 
to have nearly the same orbital characteristics as 9 comets. A further 
20 streams have been tentatively identified as of cometary origin, ' 
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reference 5> In the same reference, Whipple states that the only proven 
source of meteors is cometary, and that when more comet orbits are known 
it may be possible to link them to other streams. 

The sporadics or non-shower meteors are assumed to be dispersed 
particles from older comets which no longer travel in their original 
streams. 

The cometary- meteor model has been described by Whipple as an icy 
conglomerate and by Opik as a dust-ball. Both concepts envisage a loose 
aggregate of particles .with an overall density of <1 gm/ce, with the 
individual particles having higher densities. 

Whipple's value, for the density of -cometary meteors, associated 
with his 1957 estimate of = 25 gm/cc was 0. 05 gm/cc, reference 3- 

This value resulted in plausible luminous efficiencies of the order of 

10 J and was consistent with the frangible cometary model. He notes 
that assuming a. density of -stone, 3*5 gm/cc, led to the impossible con- 
clusion that practically all of the kinetic energy of an ordinary 'meteor 
is transformed into light. 

Other 'estimates of the density of cometary meteors obtained verbally 
from MeCrosky, Cook, and Jaechia of Harvard, ranged between O.-3 gm/cc 
to 0.6 gm/cc. 

It was agreed that the value of 0.5 gm/cc being used in the Apollo 
calculations was a good number. 

In the opinion. of Jaechia, the density of meteors associated with 
the Taurid shower is greater than that for the average cometary meteor. 
This is because the associated comet Encke is old enough to have lost 
its low density particles and the "core" particles are now predominant. 

He estimates the composition is rather like sandstone and that the' 
density is approximately 1 gm/cc-.- 

Assigning a single meteoroid density to the’ entire -magnitude range 
does not-. seem p lans ible. Whipple, in reference 3, states ‘that there is 
no evidence to justify the extrapolation .of densities' measured among 
visual meteors to the extremely faint ones. He says that the argument 
that the fragments from a large mass of low density could very well be 
of higher density is tenable. 

Whipple's latest approach to the subject .of density variation is 
given in reference 7 as the relationship, 


log io p ' - 1 * °5 " * 2li4 - ldg 10 m 


( 2 . 12 ) 
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where p is the meteoroid density in gm/ce and m is the meteor mass in 
grams. This' relationship allows p to vary from a high value of 7-8 gm/cc 
for particles of a few microns diameter to very low densities' for "bright 
meteors. 

The alternatives to this approach, adopted in this paper, are: - 

a. An all-come tary flux model with P ra 1=1 5 gm/cc or, 

b. A 'dual flux model in which 9 0 percent has a value of ' 

= 0. 5 gm/cc and 10 percent P m = 3* 5 gm/cc. 


2. '6 THE DISTRIBUTION OF METEOROID ORBITS 


Until very recently, photographic and radio evidence showed- that 
meteoroid orbits are concentrated in a narrow band within an inclination 
of 25° on either side of the ecliptic plane. 

Recent evidence obtained from the 1962 Harvard Radio Meteor obser- 
vations, up to. magnitude +10 confirms the presence of a "toroidal" belt 
of circular orbits suspected by earlier investigators. This belt, is at 
90° to the ecliptic. Up to magnitude +5, the concentration is as de- 
scribed in the first paragraph, while between +5 and +10 there is a 
mixture of both concentrations. 

The majority of sporadic and shower meteors are in direct orbit 
around the sun, that is, they approach the earth from behind while the 
toroidal meteors are retrograde. 


■5,0 HIGH VELOCITr IMPACT ON' METALLIC TARGETS 


Penetration equations for particles impacting semi-infinite metallic 
targets are of the form: 



where p is the penetration depth in the target, d is the 
characteristic dimension of -the projectile, and X^, X^, X^, X^, 'et cetera 

are physical and mechanical properties of the projectile' and the target 
material in dimensionless form. 'K is a coefficient which compensates 
for the properties which have' been ‘neglected. 
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The most reliable experimental data obtained to date has been with 
the use of the accelerated piston light gas gun, which is capable of 
1 " ' 1 " 

projecting -jr and "g glass spheres up to 22,000 feet per second. 

Various experiments have covered the velocity range from about 
5,000 feet -per second up to 22,000 feet per second using different com- 
binations of projectile and target material properties. 

The number of dependent variables involved in the penetration process 
identified so far has been limited by the necessity of varying no more 
than two or three dimensionless parameters in obtaining a curve of least 
fit to the experimental data points. The dimensionless coefficient, K, 
compensates for the parameters which have been neglected. A recent inno- 
vation in the search for a generalized penetration law has been the 
appearance of computer programs which can handle up to 27 dimensionless 
variables in a mathematical expression which will give a surface of best 
fit in a hyper-plane. The lack of reliable really high speed data points 
and the ’scarcity of information regarding the many physical and mechani- 
cal properties of the projectile and ’target materials required for the 
program will delay the appearance of a generalized ‘penetration law for 
some time. 

bn the 'meantime, the pressing needs of the Apollo program has forced 
the adoption of an interim penetration equation. 

At this point, it is necessary to say something about the current 
state of knowledge of penetration mechanics hp to the limits of present 
day projection devices. 


5.1 SEMI- INFINITE TARGETS 


Figure 10 shows a plot -of the ratio of penetration depth to the 
projectile diameter as a function of the ratio of the' impacting velocity 

and the speed of sound in the target material, V /C, , and the density 
/ P ^ 

ratio Pp/P-f- for high-speed impact into semi-infinite targets, reference lU, 


At low values of 


V /c 
P 


V 


the penetration depth varies as the 


power 


of the velocity, the theoretical predictions being -substantiated by 
several investigators. Beyond the upper limit of this region, the pene- 

2 - 

tration depth varies as the,— power of velocity. There appears to he 


a transition region in between the two slopes in which the penetration 
depth can change with velocity in almost any fashion, depending on various 
physical and mechanical properties of- the projectiles, and targets. 
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The shape of the crater and the condition of the impacting projectile 
has been determined to be a function of the impact velocity and the ratio 
of the projectile density to the target density. At low speeds, the 
crater is deep and narrow, the projectile penetrating the target without 
deformation, as the dynamic impact pressure is lower than the strength of 

lj./3 

the particle. This is the region where the V ' law applies. 


With an increase in velocity, the impact pressure is high enough 
to cause the projectile to break' up into a few large pieces. This marks 
the beginning of the transition region mentioned previously. With a 
further increase in velocity, still within the transition region, the 
impact pressure is high enough to cause an increase in the number of 
fragments and a decrease in their size. The cavity is not quite as deep 
or as narrow in this velocity range. The final stage is reached when 
the projectile and target, under very high pressures generated at impact, 
are in the fluid state and the crater shape is hemispherical. It is in 

this region that the penetration depth is proportional to 


In addition to the hemispherical crater, the rear side of the target 
sometimes exhibits a phenomena known as "spallation. " This is a circular 
fracture caused by the tensile stresses produced by the reflection of 
the shock wave from the rear face and is a feature of brittle materials. 


Reliable experimental data relating to projectile impact velocities 
exceeding a V^/c^ ratio of 1.4 on structural materials such as aluminum, 

titanium and steel, or an impact velocity of about 23,000 feet per second, 

2 

is not presently available, hence it is necessary to assume that the ^ 
power and the fluid impact phenomena are applicable at meteoroid velocitie 
It is interesting to note that tests conducted on lead targets show a 

2 i 

continuation of the — slope up to the limiting V /C. value of 4, reference 
14. . 5 - p t 


3. 2 PENETRATION FORMULAS FOR SEMI- INFIN HE TARGETS 


The penetration equations currently in use in the evaluation of 
Apollo structural reliability are developments of: 


(a) The Charters and Summers empirical formula, which is, 


R 

a 


2.28 



2/3 



2/3 


(3.1) 
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This is a modified form of the original Charters and Locke expression, 
(t) The theoretical Bjork formulas; 



p = 1. 09 . 

(m . V ' 

\ P 

> 1/3 

For an aluminum projectile 
and target 

(3.2) 

and 

p = . 606 • 

(m * V ^ 

v » py 

i l/3 

For a steel projectile and 
target 

(3.3) 


The empirical Charters and Summers expression was chosen because it 
has been substantiated by several other investigators with some slight 
variation in the coefficient and exponents, reference 11, and also be- 
cause it was determined by reliable tests conducted at the highest veloci- 
ties achieved to date, ' It should be mentioned that the light gas gun 
capable of velocities up to- 25,000 feet per second was .developed at the 
.Ames Research Center, and that Charters, Locke, and Summers conducted, 
all their tests at that facility, 

Bjork in his theoretical treatment of the hypervelocity impact 
process neglected: 

(a) The elastic waves because the stresses and particle velocities 
carried by them will always be less than those caused by the 
shock wave. 

(b) The projectile and target material strengths which are always 
greatly exceeded by the impact pressures. 

(c) The viscosity of the target material. 

(d) The conduction of heat during the period of crater formation. 

La addition, he assumed a cylindrical projectile of l/d = 1, a 
semi-infinite solid target, and the same material for the projectile 
and target, reference 12, 


The substantiating experimental evidence claimed by Bjork in the 

p/* 

same reference can apply equally well to the V law, Riney, in 

P 


reference 13, commenting on Bjork* s assumptions, reports that the de- 
pendence of the cratering process on the material strength has recently 
been demonstrated by tests carried out at the Carnegie Institute of 
Technology. He also suggests that there is evidence to show the im- 
portance of target viscosity on the impact phenomenon. 
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In the same reference, Riney rewrites Bjork's equations for aluminum 
and steel in terms, of V /C. thus: 

P t / 

\ 1/3 


Al on A l E = 2.09 
d 


'V 


V 

\l/3 


/v y 

Fe on Fe — = 1. 69 ( fr~ ) 

d V t / 


( 3 . 4 ) 


(3.5) 


in which the increased coefficient compensates for the change in form. 

Plots of p/d as a function of V /C^, using experimental data for 

aluminum on aluminum and steel on steel up to an impact velocity of 

2 

16,700 feet per second show a very good match for the — slope and a 

1 • ^ 

mismatch for Bjork's ~ slope. These plots are reproduced in figures 11 
and 12 with the addition of the Atkins and Partridge data points mentioned 
hy Bjork in reference 12. It can be seen that even at these higher 
velocities, there is no conclusive evidence- in support of the — slope. 


The newer projection techniques such as the shaped charge, exploding 
wire, the electrostatic accelerator, as well as the developed light-gas 
gun, will yield reliable test data at velocities above 33*000 feet per 
second, in the not too distant future. Until that time, it will not be 
possible to resolve the difference between the empirical and the theo- 
retical approaches to the cratering process. 


3.3 DOUBLE-SHEET THIN METALLIC TARGETS 


The current extent of knowledge of the penetration resistance of 
thin metallic double-sheet targets subjected to high-speed impact is 
primarily due to the work done by Summers at the Ames Research Center. 

Early tests showed that the rear sheet of- a double-sheet target, 
depending on its separation distance from the front sheet, could either 
be penetrated, perforated, or just dented. This led to the adoption 
of the "ballistic" limit" as a measure of the penetration resistance of 
a multi-wall target. The ballistic limit is defined as the impact 
velocity required to damage the rear sheet of a structure to the point 
where it will no longer hold a pressure differential of 1 atmosphere 
without leaks. Spallation from the rear surface of the rear sheet is 
obviously accounted for by this criterion. 
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Figure 13 is a plot of the ratio of the total sheet thickness to 
the projectile diameter as a function of the "ballistic limit, reference 
15. It is evident that there are three regions of different slopes as 
in the case of the semi-infinite slab, reference l4. 


In the first region, the total thickness of sheet required to resist 
penetration is a function of the first power of the velocity. There is 
a transition region and finally a high speed region in which the particle 
and target sheet exhibit fluid flow and the total thickness of sheet is 


proportional to the 


2 

3 


power of the velocity. 


In the low speed region up to about 10,000 feet per second, the 
projectile and the -material punched out of the front sheet forms a tight 
cluster of relatively large fragments. The major damage to the rear 
sheet is confined to a small central region with an overall scattering 
of craters. At high speed impact, above 20,000 feet per second, both 
the projectile and the target material removed from the front sheet are 
shattered into a thin shell of fine fragments. The front of the shell 
is hemispherical and travels at a slightly lower velocity than the 
original impact velocity, and after it strikes the rear sheet the material 
splashes back to the rear of the front sheet, some of’ it passing -out 
through the hole. In high speed impact on a double-wall structure with 
a sufficient distance between sheets, the hemispherical front of the 
spray spreads itself over a large area of the rear sheet. The damage is. 
more or less evenly distributed, failure generally occurring in the form 
of small cracks that radiate out of the center of the area. This indi- 
cates that the failure mode is more a pressure rupture than a perforation. 


As a matter of fact, the ballistic limit is not strongly affected . 
by sheet spacing below an impact velocity of 10,000 feet per second. 
However, it becomes _ more effective in the transition region, and above 
'20,000 feet per second the effect of spacing is much greater, reference 
13. The fragmentation spray developing a hemispherical front will 
obviously spread over a larger area of the rear sheet given a large 
enough separation distance, which, in turn, results in a lower tensile 
stress level. 


Tests indicate that the optimum meteoroid resistant double -sheet 
structure, without a filling material between the sheets, is probably 
obtained by making the front or "bumper" sheet just thick enough to 
break up the projectile and making the rear sheet the remainder of the . 
required total sheet thickness. It is also interesting to note, while 
on the subject of the "bumper" sheet, that for a given mass per unit 
area, tests indicate that the material is unimportant. 

The effect of a "filler" material such as a low-density polyurethane, 
glass wool, or any light, porous, cellular material is to increase the 
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"ballistic limit" for tie same sleet spacing and tliekness. This is. 
due to the retardation of tie particle cluster formed by the shattered 
bumper and projectile, and tests show the most effective to be the low 
density polyurethane, reference 15. 

On the other land, a honeycomb core between tie two sheets lowers 
the ballistic limit compared to a double -sleet structure laving the same 
characteristics but no core, reference 1 6 . The cells of this type of 
core prevent the hemispherical front from expanding and thus there is a 
channeling action and a consequent concentration of force on a small area 
of the rear sheet. • ; . 

For the same ballistic limit, the sum of the sheet thicknesses in 
a double -sheet target is always less than the thickness of a single 
sheet target in the high speed impact region. As a result, it is pos- 
sible to define a penetration resistance efficiency factor K as the 
ratio of the total sheet thickness in a double sheet configuration, Et f ,, _ 
to the equivalent single sheet thickness, t, or 

£t f = K • t (3.6) 


The penetration resistance efficiency factors for various double- 
wall structural combinations of sheet spacing and filler, calculated 
mainly from unpublished test, data supplied by Summers, is given in Table' 1 
As far as is known, this table has not been published previously. It 
is presented in this paper with the knowledge that the factors contained 
in it have the general concurrence of Summers. 


3. 4 APPLICABILITY OF PENETRATION EQUATIONS TO FINITE TARGETS 


Tests on semi-infinite targets subjected to high-speed impact 
indicate that as the thickness of the' target is reduced -a point is 
reached when the target is penetrated by the combined effects of crater- 
ing and spallation. The thickness of the target when this happens is 
approximately 1.5 times the crater depth, ' Thus, in designing a single 
sheet to withstand penetration at a given impact velocity, the thickness 
of the sheet would have to be at least 50 percent more than the penetra- 
tion depth calculated by. using the formulas developed for the semi- 
inf initely thick targets. Hence, it is assumed that any formula derived 
for thick targets can be used for. thin sheet provided it is modified 
to include the spallation factor. 
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Summer's equation, when modified for this target, is: 



. (3.7) 


Bjork's equations, in the form 
become: 

Al on Alt t 
Fe on Fe: t 


given by equations 

(3*2) and (3*3), 

=» i. 6k ( ra p v p 

(3*8) 

= 0.91 (mV Vfr 

(3* 9) 

\pp / 


Equations (3-8) and (3*9) ma y be approximated by a, 'single equation, which 
is assumed applicable to all combinations -of target and -projectile 
materials. The generalized equation is: 


P p l/? •, P t • (Vp ) 1/3 (3 - 10) 


A further modification introduces the target material sonic' velocity 
into the equation, Which becomes: 


5.83 


1/3 


-1 



(3. 11) 


This equation is valid for all materials in which the sonic velocity is 

8 

5 km/sec, or for which the Young ' s modulus to density ratio is 10 . The 
materials falling in this category are the usual structural materials, 
aluminum, titanium, magnesium, and steel* 


would be: 


For beryllium, (c^ = 12.6 km/sec and - 1,835 gm/cc j , the - equation 
be: ' ' ■ 

, / XT \1 h 

1/3 


t = 7*95 p. 


-l 


V 

m . -£ 

P C 4 


= h.Jh 




( 3 * 12 ) 
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; 5.5, the PENETRATING- METEOROID FLUX - 


As suming a spherical meteoroid. Summer's penetration equation for 
finite sheet, equation (5*7-Xf can .he. rewritten as: • 

t - k.zk p h3 . m 1/3 . ( Ji_ 
m I •• Y e P t 

where 

t = The equivalent single sheet thickness of a structure 
just capable of holding l‘psi pressure- differential, 
cm. 



V - Meteoroid impact velocity, km/sec 


m 


Cj. = Sonic velocity of structural material, km/sec 


= Density of meteoroid, gm/ec 


Pt 


= Density of structural material, gm/cc 


m = Penetrating meteoroid-mass, gm. 
3? * 


From equation (5. 15 )* 


m = 


El. 


P 16 . 2 b p . p 

m t 


—2 


V /C 
m t 


( 5 - 1^0 


The general form of the meteoroid flux-mass relationship- is given - 
in equation (2. 11 ) as, 


N 


a < 



a, and 
and ST is 
stituting 
rate, . 


0 are constants, m is the meteoroid mass in grams, 

the. flux rate in gram impacts- per- unit area per day. Sub- 

m for m results in an expression -for the penetrating flux 




2k 


A similar expression can “be -obtained by using the modified Bjork 
equation, ( 3 . 11 ), that is. 



(3.-16 ) ' 


Thus, for every structural component exposed to the meteoroid environ- 
ment, there is a penetrating flux-rate, N , which is the number of 

penetrating meteoroids per unit area and time of mass and heavier .- 

Figure it is a plot of BT^ as a function of t for an aluminum 

structure, combining Summers and Bjork' s equations with the current 
Apollo -meteoroid flux-mass expressions. 


k. 0 PROBABILITY OF PENETRATION 


Having determined the penetrating meteoroid flux for a given 
equivalent skin thickness and material, the total meteoroid population 
can he considered as consisting of only two types — penetrating and 
nonpenetrating. ' - - 

If a trial can have one of two mutually exclusive results, 
(penetration and nonpenetration) and if the probability p of a 
penetration is constant over a series of n independent trials, then 
the probability of obtaining r penetrations is, x 


P 

r 


nl 

rj (n-r ) ! 


P 


n-r 

1 


(^. 1 ) 


where, q = (l-p). ; ’ 

The probability distribution of r is called the Binomial distribution 

since the probabilities can be obtained by expanding (q+p) n by the 
Binomial formula. 


In the case of a spacecraft in the meteoroid environment, a trial 
is considered an impact and the probability of a penetration in a single 
impact is: ^ 


( 4 . 2 ) 
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where N. is the total meteoroid flux and U is' the penetrating flux 
t - p 

for the spacecraft. 

The -value' of p is assumed constant during a mission, as the effect 
of removing the meteoroids that strike a space vehicle- from. the total 
flux is insignificant. As an example, say the flux -mass equation is 
-9 . 

N Xm = lXlO gm impacts per square foot per day. • Then- for a 
x« 

10,000 square foot sampling area, the number of impacts up to and in- 

-12 ' 

eluding a mass of 10 grams is, 

N = 1X10 7 per day 
"G 

For a given equivalent skin thickness of a space vehicle, the number 
of penetrating meteoroids for the same sampling area is,, 

KT = 10 per day- ■ 

P 

Then at the first trial, p. = 1X10 ^ 

After 10° trials, all of which were removed from the environment, 
and assuming one was a penetration, 

jp = — 2__ = ^ 10 -6 

9X10 6 


^ • ' 

If* the 10 trials had included *1- penetrating meteoroids, then. 


p a = 6 . 67X10 -7 

9X10 

li is always very much smaller than N, , hence the assumption- of 
p t 

a constant p is valid. . ■ 

The "mean" or the "expected -number of penetrations" of the proba- 
bility distribution of ■ r‘ is the product -of the number of trials, h, 
and the probability of a single penetration in a single trial or impact, 
p, that is. 

Expected 1 number of penetrations, r g = n-p. 


( 4 . 3 ) 



(4.4) 


The number of impacts , n - * A • t . p| 

where N. = total meteoroid flux, gram impacts per unit area per day. 

W * 4 

A = vehicle surface area 

T = time in the environment in days 

T\ s earth shielding factor, which varies between 0.5 and 1.0 
The substitution of equations (4.2) and (4.4) in (4.3) results in, 

r = n-p = N • A* r- T) (4.5) 

e p 

When the probability of a single penetration in a single impact is 
a "rare event", (p is very small) and the number of impacts is very 
large, the use of the binomial probability distribution relation in its 
exact form becomes unwieldy. Poisson’s approximation to the binomial 
distribution, 

p =Ae£ . 6 - a -r (4.6) 

r r! 

Substituting for n. p from equation (4.5) gives, 

P . (v A "f 

r rl 

where r is the number of penetrations. 

The probability of no penetration is , 

P a e" n * p 
o 

or 

loS 10 P o = “°' ^ n * p = 

This equation can be used either as a means' of designing to a given 
probability of zero penetration, or to evaluate the value of P for 
a given structure. 

.The probability that the number of penetrations will not exceed 
a number m is the summation of the probabilities of each occurrence ' 


■ -n .A.-r.'n 

e p 11 


( 4 . 7 ) 


-0.434 IT .A. T.p 

*n ■ 


(4.8) 
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up to m, including zero, that is 


P 

r < ra 


=. P + P n 
o 1 


+ P 


2' 


+ P 


r=m 



or 


**N , A. t * Ti 

P = e p X 

r £ m 


r=m 


1 


N . A. t . 1 \ 
P 


r! 


(4. 10) 


The probabilities of exactly 0, 1, 2, 3’, and 4 'penetrations and of 
not exceeding 1, -2, 3, or 4 penetrations for n.p values ranging from 
0 to 1 is given in table 2. 

The binomial- distribution,- when presented • in the graphic- -form, -is 
represented by a histogram and not a continuous curve. ■ Events -are • 
measured along the abscissa, probabilities along the ordinate,- and' the- 
histogram. is constructed so that each rectangle is- of unit width and 
centered about r = 0, 1, 2, 3 , and so- on. It follows, that the prob- 
ability of each event r is also the area of the rectangle on r, 
hence the sum of the probabilities up to any value of r is the>'area 
of the histogram up to that, ordinate. The histogram- for' -an expected 
■number or "mean" of 5 is shown in figure 15* 

The total area of a probability histogram can never exceed unity 
because by definition the summation of all the admissible values of 
r, including zero, of the "probability f unct-ion,, " equation -( ; 4 .6--), is, 

r=s 

V e -nP = iaL]I = i (4.ii) 

r=0 


3. 0 RELIABILITY AMD COUFIDMCE LEVEL 


Although the probability of no penetration is generally used as the 
criterion for the structural integrity of a space vehicle exposed to the 
meteoroid environment, it is not a -measure of reliability. The reliabil- 
ity of a component is usually determined from -a statistical sample of 
identical tests, which is obviously not a practical proposition in the 
case of a space vehicle, 
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A method of calculating the reliability , together with an associated 
confidence level based on the expected number of penetrations over a 
very large number of identical missions, will be described in this section. 

For a situation where n.p is the most probable number of penetra- 
tions, it follows that, 

n.p = the most probable or average failure rate 
= F 

From this the most probable reliability is, 

R - 1 - F (5*1) 

It is evident from figure 15 that in dealing with the most probable 
number there is a good probab ility_of exceeding it. Thus the failure 
rate may be equal or greater than F approximately 50 percent of the time, 
and the actual reliability may be equal or less than R approximately 
50 percent of the time. The confidence level, C, in the reliability 
being R or more is therefore about 50 percent. 

» 

In order to achieve a reliability R with an associated high 
confidence level the system must have a most probable failure rate less 
than (l-R), that is 


F < F = (l-R) (5.2) 

The relationship between F and F may be illustrated by the 
two probability distributions shown in figure 1 6, which differ only in 
the number of trials (mission time). 

To attain a failure rate in which one has 97 • 5 percent confidence, 

1 2 

the average failure rate is approximately ~ F when np = 5 and — F 
when np =» 15 * 

This .ratio will be called k, _ Thus, 

F = k.F 
a k. (np ) 

and R - 1 - F 

= [l - k, (np)] 


(5-5) 
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Values of k for various confidence levels are plotted in figure 
17 as a function of (a. np), -where "a" is a large number of missions such 
that (a, np )» 1 . 

It may be seen that for C a 0, 5> k-»l. 0 as may be expected. 

Also, for, 

C = 0.90 k-»1.3 

= 0,95 k-»l.4 

=* 0.975 k-*1.5 

Por reliability purposes, the generally accepted confidence level 
C = 0.975 will be adopted, in which case k a 1 . 5 , 


Thus , ' ‘ r n . 

R * [l - (5-M 

when C = 0.97-5. 


This means that, ’in order to quote a reliability in which one has 
a high confidence level, this reliability must be based on a. failure 
rate which is 1 , 5 times the average failure rate. , Conversely, the average 


failure rate must be ■=■ of the guaranteed failure rate. 

5 


6 . 0 RELIABIL3TY - EQUIVALENT SKIN THICKNESS RELATIONSHIPS 


The principles involving the evaluation of the shielding require- 
ments for a vehicle in the meteoroid environment have now heen established, 
All that remains is to directly relate reliability to the equivalent 
skin thickness in terms of flux parameters, penetration equation, and the 
number of failures. 

Substituting equation (4.5) into', equat ions (3*15) and ( 3 , 16 ) gives: 
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(Id) Using Bjork's equation, 

np = a. A. t« T] 

Thus, by equation (5.l) 


(t) 5 , P t A 

200 Rn 


■ V , 

IllJ 


(l - k. a. A, t.T] 

A 

l 

76.24 p' .v 

“ m m J 


and 


using Summers' equation 

IP 


R = /l - k.a. A. t.T] 


<■ 

__ j 


m . m 
using Bjork’s equation 


( 6 . 2 ) 


(6.3) 


(6.4) 


The penetration equations considered are examples of the most generally 
accepted empirical and theoretical expressions of the cratering phenomenon. 
Other forms of penetration equations are equally amenable to the procedure 
established in this paper for the ultimate relationship of equivalent 
thickness, T to the protective reliability, R. 


A. 1 MSG DESIGN EMVJROMMENTS 


The design environment currently in use within MSC divides the total 
flux into two components in the following manner: 

a. Asteroidal: ■ Flux =» 10 percent of the total 

m =1 gram 
o 

P m * 3.5 

3 = -l.oo 

^10 2 
a ** 4.5 X 10“ gm impacts/meter day 

V - = 30 km/sec, or as in figure 9 
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m 

0 

Pm 

3 

a 

V 

m 


Thus, by equation (2. 11 ) , 

t i J T 

, * ' *<] 0 ' - 2 

If. m = 4. 5 X 10 “ gm impacts /meter day 

Or in engineering units , ... : 

N. m = 4.2 X 10- gm impacts /feet day 

b. Cometary: Flux = 90 percent of the total’ 

= 2.0 gram 

= J 0, 5 gm/cc 
= -1 

a 8.2 X 10 ^ gm impacts /meter^day > 

= 30-ton/sec, or as in figure 9* 

U.m « 8. 2 X 10 “ gm impacts /meter day 
= 7- 66 -X 10 gm impacts /ft^day * 

The design environment is equivalent to a modified form of Whipple’s 
1957 flux -mass relationship, in which the value of m =2.5 grams, as 

suggested in reference 7, and a meteoroid density of 0.5 .gm/cc is used 
instead of 0.05 gm/cc. Thus the modified 1957 Whipple flux mass re- 
lationship is: 

c. N* m = 1. 134 X 10“ gm impacts /meter aay 

a 1.06 X 10 gm impacts /ft^day 

and 

"^m ~ 50 km/sec, or a's- shown in figure ■ 

The equivalence of (c) to (a). + .(b). .can be -demonstrated as follows: 
From equation. (3 . 15) with p = -1 


and 


If = a . 
P 


76.24. p .v- 

* r m m 

(t) 3 .p t 2 .o/ 


, where t is in cms. 
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Substituting for = 30 km/sec, a 5 km/see and p t = 2.78 gm/cc. 


35^*24 a, p. 


N = 
P 


m 


(t) 3 


(1) 


n . 1T n 

For (c), substituting for a = 1.13^ X 10 ~ , t * 10“ cm, and 


P =0.5 
r m 


cc 


N p = 2. 01 X 10-5 p ene t rat ions / f t^day - 

_ “1 

Also, for (a) + (b), again using t s 10 cm, and the appropriate 

values of a and p . 

r nr 

- 3 - 5 ^ X v? (^.5 X 10“ 10 x 3,4 -t- 8,2 x 10" 9 x 0.5) 

= 1. 99 X 10 "5 penetrations/ft ^day 

Hence, in terms of penetrating flux, and also np, the dual flux 
may be replaced by the 100 percent cometary model, (c). 


A. 2 CALCULATION OF P MD RELIABILITY, R 

o 7 

The basic data for a space -vehicle component on a lunar mission is 
as follows: 

p 

Surface area, A a 520 ft 

Lays in meteoroid environment; t a ik 

Construction: An aluminum, , double -wall without a filler, with. 

equal sheet thicknesses of .017 in., 1.0 in; apart. 

•Design environment: N.m a 1. 06 X 10 ~ 9 gm impacts /ft^day 

P m a 0. 5 gm/cc 

V a .30 km/see 
m ' 
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From equation (3*^)’> 


equivalent thickness , t = 


2t f 

IF 


Using table I, K = 0,37 for a 1.0 in. double -wall, no core 

: t . ’ - ' * ’ - 

*• .37 

= 0.092 in. 

On substituting in equation (6. l), for a = 1.06 X 10 ^ 0 = “!#p m = 0.5 ^m/cc, 

p ~ 2.J8 gm/cc, V = 30 Kin/sec, and C = 5 Kin/sec, and using an earth 
shielding factor of 1.0, the 

—7 A. t 

expected number np = 1. 88 X 10 x rv , where t is in cms (2j . . 

■ " . • • (W ■ • ■••• • ■ 


Converting to inches, , 

np = 1.1^7 X 10- x^r* 
■ ■ (t) 3 


(3) 


Using equation (6,2), a similar expression for the expected number 
is obtained, that is, 

np. a 2. 96 X 10*"® x ■— ’ -■-• j , - for t in cms. 

( t ) 3 

and _ q . 

np » 1.8l X 10“ y x " for t in inches (U) 

(t ) 3 

Figures 18 and 19 are plots of np as a function of t for various 
values of the area-time product A', t for both Summers 1 and Bjork’s 
penetration 'theories. 


Hence, when t = 

.092 in. 

. 0 . - 

and A. T « 7 >290 ft. days. 

a. 

Summers: 

£ 

11 

k 

-u 

b. 

Bjork: 

np = .017 


The probability of no penetration P is obtained from the curve 
given in figures 21 and 22. 0 



Thus for the Summers' case, 

P = o, 899 
o 

and for Bjork's, P q = 0 . 983 

The corresponding reliabilities and confidence levels obtained as 
. described in section 5 are as follows: 



np 

P 

0 

C 

R 

(a) Summers 

.107 

.899 

50.0^ 

.893 

- 



97. 5^ 

.83: 9 

(b Y Bjork 

. 017 

.983 

50 . 0 ^ 

.983 




91.% 

‘ ' * 975, 

If the void between the sheets 

were filled with a 

honey-comb 



K a 1.00 



and- 


t a 0*034 

in. 



np . ^ 

P 

0 

C 

R 

(a ) Summers 

2*12 

.120 

50*5^ 

0 

' 

1 


97 .% 

0 

(b) Bjork 

•33^ 

.716 

50.% 

. 666. 



■ 

91.% 

\499 


' A, 3 DESIGM3 TO A GIVM P’ ORB 

• ■ . o 


A* 3.1. Based on P q « 0*99 — 


n 

Prom equation (4.8), log^g P q -» - ‘ 1 np 
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The equivalent single sheet thickness in aluminum is obtained’ from. 

O 

figure 18, using the Summers curve, and A.t = 1,290 ft days. 

Thus, t = 0,202 in. 

. Using table 1, and = K. t, and considering a 2. 0 ,f double-sheet 

construction without a core, 

2t f = .23 X' .202 
= 0,0465 in. 

The sheets could he equally divided into two 0, 024 in'. or the' front 
or bumper sheet may be made only as thick as the diameter, of the pene- 
trating meteoroid, and the rest' of the material added to, the' rear sheet. 

In order to calculate the diameter, it is ' necessary to obta'in the 
penetrating mass from equation (3? l4), in inches. 

_ 3 p 

M ■ H 

m r 

p -4, 65 p . (v /c . r 

m m t 

.The diameter is- calculated hy. assuming' a spherical meteoroid and 
the use of the appropriate density. 

For this structure., R = 0. 99 for C = 50^ ‘ 

and R «. 0. 985 for C = 97. % 

A. 5. 2. Based on R g >0. 99 and’ ~C = 97. 5$. — 

From equation (5,5) R = jjL - k(n.p,)J 
when ’ C = 97*5$, k » 1.5 

0.99 « [l - 1.5 n *pJ 

n. p =0. OO67 


and 
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As before, the value of t can be obtained from figure 17, using the 
same criteria. 

t = 0.232 in. 

and 

rt f « .23 X .-232 

= .033 in. 

Designing to a confidence level of 0. 975 and a reliability of 
O.99 results in a l4 percent increase in weight. 

However, the increased confidence level in the reliability of the 
structure is sufficient compensation. 

The question of material selection is covered in section B2. 'Sheet 
thickness and weight ratios relative to aluminum are given therein for 
the other commonly used materials. 

As an example, the percentage increase in weight due to the increased 
confidence level can be offset by using magnesium instead of aluminum, 
providing its use is permissible, 

A plot of np as a function of reliability for confidence levels 
of 97-5^, 95^, and 90$ is shown in figure 20. 


A, 3.3. Designing to a. controlled number of penetrations — 

Sometimes it becomes necessary, due to weight limitations, to permit 
no more than a given number of penetrations within a specified time 
period. The procedure to be followed in this case is identical to that 
outlined in the previous sub-section, A3. 1. 

Figure 23 is a plot of .np ■ as a function of the number of events 
that should not be exceeded, for probabilities of 0,99, 0. 995, and 0.999* 
Thus, if a certain design calls for a 0.99 probability of not exceeding 
two penetrations in 100 days, then from the curves , 

np a 0. 4-5 


Using the Summers' 


criterion, equation (2), 



1.147 X 10 


-8 


x A. T 


np 


9 


where T = 100 days. 
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In this case, the concept of "reliability". is -no- longer valid. A 
structure that is designed to be penetrated is not "reliable" in the 
sense previously outlined for no penetration. 

There is no need to specify anything other than the confidence that 
the probability of not exceeding n penetration in t days is 0. 99 or 
0.999* 

The confidence is obtained from the lower curve of figures 21 and 
22 which, for the example being considered, is 0.537* 

Design charts for the dirept determination of t for P ^ =0.99 

and Summers * penetration equation are given in figures 2 b and 25 • 

In one case the number of penetrations is limited for the' total 
mission time t and in the other for various lesser time periods. 

, 2 

As an example, a space-vehicle of lo , 000 ft is expected to remain 
in the environment for- 500 days. 

6 2 

Then total At = 8 x 10 ft days' 

Based on the philosophy of no more than one penetration every 
50 days,' t from figure 25 is 0. 5^5 in. 

If an average of, no more than one penetration in 30 days were based 
on an actual six in 180 days, then t a 0.24 in. 

Should the average of one in 30 be based on an actual 17 throughout 
the total mission time of '500 days, then t is 0.175 in. approximately 

Obviously, the actual design is based on the' requirements of the 
vehicle or component to be protected. 

A thickness of 0. 3^2 in. insures that there will be no more than 
one in 30 days, whereas the value of = O.-I.75 in. being based on the 
average does not preclude the possibility that all 17 penetrations will 
occur in the last 30 days of the mission. The latter condition, being 
much more hazardous to both' -crew and vital sys.tem components, could not 
be tolerated. 



B. 1 MATERIAL COMPARISON 
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B, 1, 1 Sonic Velocity in Structural Materials — 

0 .(iif 

where - 2 

E a Young's Modulus, lb/in 

p a Density, Ib/ir P 

U 

g = 386 , ■U- in/see^ 


Material 

\ 

\ x 10 ' ~ 6 

(ft/sec) 

Cj. (km/sec ) 

Magnesium 

0. 065 

6*3 

16,350 

5 

Aluminum 

0.1 

10 

16,350 

5 

Titanium 

0.16 

16 

16,350 

5 

Steel 

0.30 

30 

16,350 

5 

Fiberglass 

O.O 63 

3 

11,300 

3*46 

Plexiglass 

0.043 

0.43 

5*170 

1.58 


B.2 SINGLE THICKNESS /TO WEIGHT RELATIVE TO ALUMINUM 


B. 2.1 Summers' Equation — 

From equation (3*13)* converted to inches, 


t = 1 . 67 p, 


„ i/5 


2/5 


in 


m 


or 


• 16 (Vt) 


m ■ -p 
-2/3 


* \C*j_ * 


V p t. 


cms 


(5) 
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Material 

t.Cto). 

t ^ t Al 

,P t (it /in 3 ) 

w (lb /in 2 ) 


v/wai 

Mg 

9.6 x io“ 3 k 

1. 33 

0, 063 

6 . 25 x- 10“^ 

K. 

0.87 

A 1 

7. 2 x 10 -3 K 

1.00 

0. 100 * 

7. 20 X 10 

K 

.1. 00 

Ti 

5 . 28 x 10" 3 K 

0 . 73 

0.160 

8 . 45 x I 0 ~ k 

K 

'1. 17 

St 

3.55 X 10“ 3 k 

0.49 

0. 300 

_k 

10,65 x 10 

K 

1. 43 

KG. ■ 

12 -. 5 x 10 ** 3 K . 

1.73 

•0.063 

-4 

7. 87 x 10 

K 

1. 09 

P.G. 

27. 1 x 10“ 3 K 

3.76 

0.043 

-4 

11.65 X 10 

K 

1. 62 


B.2.2 Bjork’s Equation 


From, equation (3.11), in inches 


t = 2.30 K.p t " 1 .C t -“ 1 ^ 5 

' (6) 

for Al, .Mg, St, and Ti 



Material 

t (in) 


^Al 

P t (lb /in 3 ) 

w (lb /in 2 ) 

W/W A1 

Mg 

13. 9^ X 10" 1 

K 

' 1. 54 

0. 065 ■ 

. 9.07 X 10 " 2 -K 

•1.00 

A1 

9. 06 x 10” 1 

K 

1.00 

0.100 

9. 06 X 10 “ 2 K 

1.00 

Ti 

5 .66 X 10 ** 1 

K 

0.625 

0. 160 

9.06 X 10 ~ 2 K 

1. 00 

St 

3. 02 x 10 _1 

K 

0 . 333 

0.300 1 

9. 06 -X 10~ 2 K 

1. 00 
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r 

0 

1 

2 

3 

4 

n.p 

p o 

P ! 

P 2 

P_ 

3 

1 

P 4 

0 

1.0 

- 

- 

- 

- 

.01 

. 99005 

.00990 


- 

- 

.10 

* .90^4 

.09048 


.00015 

.00004 

.20 

.81873 ■ 

.16375 

.01638 

.00109 

.00006 

'.50 

.60653 

.30327 

. 07582 

.01264 

.00158 

.80 

.44933 

■35946 

.14379 

.03834 

.00767 

1.00 

. 36788 

f 

.36788 

.18394 

. 06131 

-01533 

1.10 

; .33287 

.36616 

.20139 

.07384 

.02031 


TABLE II.- PROBABILITIES OE 


<1 

<2 

<3 

-=d* 
t V 

P + p 

0 1 

V P l + P 2 

p + p + p + p 

0 1 2 3 

p + p + p+ p + p, 
0 1 2 3 4 

.99995 

1.00000 

1.00000 

1.00000 

.99532 

. 99984 

•99999 

1.00000 

.98248 

.99886 

.99995 

1.00000 

.90980 

. 98562 

. 99826 

. 99984 

.80879 

. 95258 

. 99092 

.99859 

.73576 

.91970 

. 98101 

. 99634 

. 69903 

. 90042 

. 97426 

• .99457 


AND EOT EXCEEDING r EVENTS. 
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BALLISTIC LIMIT, FT/SEC x 10 “ 3 


Figure 1J* - Penetration parameter vs velocity for finite 

double -sheet targets. 
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